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Investigation of Models for Large-Scale Meteorolofietl 
Prediction Experiaents 

During the past year the City College group continued its collaboration 
with the Goddard Institute for Space Studies (GISS) in the investigation 
of climate simulations generated by the GISS coarse-nesh climate model 
(Hansen et al., 1980). 

Graduate students employed on the project were Zaphiris D. Christidis 
(through May 1980) f, Michael Dennis (through June 1980), Charles Cohen and 
Peter Wu. 

Two technical reports have been distributed ("Spherical harmonic 
analysis of a synoptic climatology generated with a global general circula- 
tion model" by Z. D. Christidis, and "Eigenvector analysis of some observed 
and model-generated climatological fields" by M. Dennis and J. Spar), A 
paper by J. Spar ("Prediction experiments with a coarse-mesh global model") 
was published in the Proceedings of the Fourth Annual Climate Diagnostics 
Wv^fkshop, Madison, Wis., October 16-18, 1979 , NOAA, pp. 413-423, March 1980, 
and a paper by Z. D. Christidis and J. Spar ("Spherical harmonic analysis 
of a model-generated climatology") was accepted for publication in the 
Monthly Weather Review (February 1981). 

The remainder of this report consists of a paper on a series of climate 
simulation experiments conducted with the coarse-mesh climate model by Spar, 
Cohen, and Wu for the purpose of better understanding the role of initial 
and surface boundary conditions on the model -generated climates. At this 
time only preliminary results of these experiments are available, and 
statistical evaluations have not yet been carried out. The experiments 
themselves are still in progress, and will be completed within the next 
few months. 


Cliaate SiBulation Experiments with a Coarse-Mesh General 
Circulation Model ; Prelislnaiy Results 


J, Spar, C. Cohen and P. Nu 


Introduction 

General circulation models (GCM’s), including the GISS coarse- 
nesh (8° x 10°) climate :nodel (Hansen, et al., 1980), are often used 
to perform "prescribed change" experiments in which some alteration 
in solar radiation, surface boundary conditions, or atmospheric com- 
position is specified and the atmospheric response, i.e. the effect of 
the prescribed change on the weather or climate, is calculated. In such 
experiments, the specified change is usually some small perturbation of 
the basic constraints on the system, and the computed response of the 
model may be weaker than the background noise. The analysis of these 
experiments would undoubtedly be aided by a better understanding of the 
ways in which the primary climatological controls combine to generate 
the basic climatic state. In this study, an attempt is made to compute 
the contributions of various surface boundary conditions to the monthly 
mean states generated by the 7-layer, 8 x 10 GISS climate model (Hansen 
et al., 1980), and also to examine the influence of initial conditions 
on the model climate simulations. 

For the purposes of the study, such obvious climatic controls as 
the shape and rotation of the earth, the solar radiation, and the dry 
composition of the atmosphere are fixed, and only the surface boundary 
conditions are altered in the various climate simulations. The model 
(version 660) is operated at a fixed solar declination, specifically 
that of January 15, initialized on January 1, and allowed to simulftte 
15 or more successive Januaries without going through the annual cycle. 
Monthly means are computed for each January and for ensembles of all but 
two (or more) Januaries, the first two (or more) months being discarded 
as transients. 


The Bodel is initially forced by the sea-surface tea^erature (SST) 
field alone (plus, of course, rotation, etc.) by eliminating all land 
and specifying a tonally uniform SST pattern with the value at each 
grid latitude equal to the zonal average of the climatological monthly 
mean SST, but with sea ice caps (3m thick) over the Arctic and Ant- 
arctic. For this "water planet" experiment (001), the model was 
first initialized with monthly mean zonal values of the atmospheric 
history variables, including winds, taken from an earlier 5-year model 
climate simulation (Chiistidis and Spar, 1980), exeunt for the surface 
pressure which was assigned a globally uniform initial value of 1010 mb. 
This initialization was adopted for the purpose of reducing the "spin up" 
time in the experiment on the assumption that the model would soon "forget" 
its initial conditions and adjust rapidly to the boundary conditions. How- 
ever, early results of the experiment indicated an apparent tendency for 
persistence of the initial state (except for a redistribution of the 
surface pressures and some alterations of the winds), probably as a con- 
sequence of the constancy of the forcing, the absence of an annual cycle, 
and the zonal symmetry of the surface boundary conditions. Indeed, as 
shown below, there was enough similarity between the initial state and 
the model-generated climate in this perpetual January "water planet" 
experiment that the possibility of initial conditions having a dominating 
influence on the model-generated climate could not be ignored. It was 
therefore deemed necessary, after all, to conduct a "water planet spin up" 
experiment (000) in which the model is initialized with globally uniform 
mean values of sigma-level specific humidities and ten^eratures, as well 
as a constant surface pressure and zero winds. 

Following the two water planet experiments, flat continents are 
placed on the earth, with zero elevation above sea level, zero water 
storage capacity, and uniform surface albedo, in order to assess the 
thermal influence of the land, while the zonal distribution of sea ice 
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and SST remain unchanged over the oceans. The model is again ini- 
tialized with the same state of rest as in the water planet spin up 
experiment, and is again run for a period of at least 15 Januaries, 
the first 2 (or more) months being discarded. The initial surface 
albedo of Greenland and Antarctica is that of glacial ice, but snow 
is allowed to alter both the sea ice and continental albedos as com- 
puted by the model. The climatic effects of the flat dry continents 
as heat sources and sinks may then be evaluated by comparing the 
’•flat continent" model climatology (002), confuted by averaging all 
but the first 2 or more Januaries, with that of the water planet (000). 

The effect of mountains on the global model climatology is com- 
puted by repeating the last run, but with the correct model terrain 
elevation restored at each grid point. However, in the "mountain" 
experiment (003) the model is initialized not only with a state of 
rest, but also with a completely dry and isothermal atmosphere, so that 
it is required to generate its own humidity and temperature distribution. 
The ground albedo is still kept uniform over the continents, except 
where snow is calculated, and the water storage capacity is zero in this 
mountain run. 

The climatic effects of variable surface albedo and water storage 
capacity over the continents are investigated through a second perpetual 
January mountain experiment (004) with the model, while the influence 
of departures from zonal symmetry of the Sfff field is examined by 
comparing the results of the previous run with the perpetual January 
climatology generated by the complete model using the actual January 
climatological SST's. 

This report presents some preliminary results of the experiments 
conducted thus far. A more complete analysis of experimental results, 
including statistical evaluation, will be reported later. 
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Bcperinents with a zonally symmetric water pl/^iiiet 

In the first of the two water planet experiments (001) the 
coarse mesh 7>layer model was initialized with zonal mean values of 
temperature, humidity, and wind for the month of January derived 
from an earlier model climate simulation based on 5 annual cycles, 
but with a constant surface pressure (1010 mb). During the first 
month of the simulation, a rapid mass redistribution occurs, undoubt- 
edly in response to the initial wind field. The 13-year mean meridi- 
onal distribution of surface pressure (P^, in mb) shown in Fig. 1 
evolves quickly, with a near-equatorial low, subtropical or mid- 
latitude highs, subpolar lows, and an Antarctic high. At the same time, 
the zonal wind distribution undergoes a moderate alteration, as illus- 
trated in Fig. 2, which shows the meridional profiles of the mean 
zonal wind ( 0 , in ms"^), integrated over the 7 layers of the model, 
both in the initial state (dashed) and as generated by the model and 
averaged over the last 13 Januaries (solid). Within the first 2 months 
of the simulation, the tropical easterlies, which initially extend from 
8N to 23S, shrink to the zone 8N to 8S and diminish in speed, 
while the strongest westerlies, initially located at 47N and 39S, 
shift closer to the Equator and stabilize as shown in Fig. 2. 

The meridional wind component (V) also changes as shown in Fig. 3 
illustrating the meridional profiles of V in the lowest layer of the 
model (Vj^, in cm S'^) both as given in the initial conditions (dashed) 
and as generated by the model and averaged over the last 13 Januaries 
(solid) . Initially, maximum low-level northerly (negative) winds are 
centered near 20N, while the strongest southerly (positive) winds ap- 
pear at 31S, with a secondary southerly maximum at the Equator and 
an apparent double intertropic^il convergence (ITC) . Within one month 
the northerly maximum shifts equatorward and the principal southerly 
maximum moves to the Equator and amplifies, producing one dominant ITC 
at 4N as shown in Fig. 3. The resulting distribution of vertical 
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velocity (W, in wn S”^), integrated over the 7 layers of the sodel 
and averaged over the last 13 Januaries, is shown in Fig. 4, where 
the Baxinum upward (positive) notion at 4N is the dominant feature 
and the cause of a precipitation maximum. Also noted in Figs. 3 and 
4 is the fact that, although latitude 20N is initially the location 
of the belt of strongest northerly winds (the northeast trades), and 
Hence almost zero vertical motion, the model generates maximum diver- 
gence, and hence maximum subsidence, at that latitude, completing a 
Hadley circulation that develops in the first month. Another radical 
alteration of is seen in Fig. 3 in the Arctic, where an initially 
southerly wind at 70N reverses (again in the first month) and re- 
mains northerly in the 13 January average. 

Deviations from initial conditions are somewhat less dramatic in 
the vertical distributions of wind and temperature. In Fig. 5, showing 
vertical profiles of the initial and 13 year average model-generated 
zonal winds (U) at latitudes 39N and 39S, plotted against pressure 
(P) on a logarithmic scale, the model exhibits a marked decrease in 
the westerlies. The lapse rate of temperature (not shown), on the other 
hand, exhibits very little change from the initial state at all latitudes, 
except over the open ocean south of the Arctic ice cap where the sur- 
face air temperatures rises quickly from an initial value of - 25C to 
the water temperature of - 1C. 

As expected from the zonal symmetry of both the initial atmospheric 
state and the surface boundary conditions, the horizontal fields of the 
model-generated climatic variables display mainly a zonally symmetric 
character when averaged over 13 Januaries, although waves and cells are 
found on the individual January maps. For exan5)le, the third mean 
January sea- level pressure field con^juted by the model (Fig. 6) ex- 
hibits two or three distinct Arctic lows and three or four subtropical 
high cells in the Northern Hemisphere. However, the 13-month mean 
Jaunary field of sea- level pressure in Fig. 7 is predominantly zonal, 
as are those of the layer temperatures for 850-1000 mb and 700-850 mb 
(Figs. 8 and 9) and the geopotential heights of the 700 mb and 500 mb 
surfaces (Figs. 10 and 11). The last four fields, in fact, remain 
nearly constant during the 15 months of the perpetual January simu- 
lation. 
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Mean neridional cross-sections of temperature (T) and zonal 
wind (U) for the 15th month of the first January water planet ex- 
periment, shown in Figs. 12 and 13, exhibit some of the defects of 
the model simulation. The tropical tropopause is too low and too warm, 
the Arctic tropopause is too high and too cold, and the subtropical 
jet in the Northern Hemisphere is much too high. However, the tropo- 
spheric structure of the model atmosphere is not unrealistic. 

While there are clearly sigt?ificant departures of the model- 
generated winds and pressures from the given initial conditions, there 
also appears to be considerable persistence of the basic atmospheric 
structure and circulation* albeit in a modified form. Furthermore, 
it appears likely that the alterations of winds and surface pressures 
computed by the model may be due largely to the initial inconsistency 
between the given wind field and the constant surface pressure, as 
indicated by the rapid mutual adjustment of these two fields. Thus, 
there is still a question as to what the model- generated climate would 
be if the simulation were initialized with a different set of initial 
conditions, or even allowed to "spin up" from a state of rest. To 
answer this question, the water planet model was reinitialized with a 
set of constant ten^jeratures and specific humidities at each level 
(specifically, the global averages of the initial 5-year model clima- 
tology), the same constant surface pressure as before, and zero winds. 
Again the model was run for 15 months in the perpetual January mode, 
and the output was averaged over the last 13 Januaries. 

The initial adjustment of the model atmosphere in the "water 
planet spin up" experiments (000) can be described briefly in terms 
of the mean state for the first month. Surface air temperatures quickly 
approach the SST's over the open ocean, while in the Arctic both the 
sea ice and surface air cool to about - 10° C, At all tropospheric 
levels a realistic meridional temperature gradient begins to emerge, 
while in the stratosphere temperatures increase monotonically from 





North Pole to South Pole. A realistic neridional gradient of specific 
huoidity begins to appear in the lower troposphere, with very high re- 
lative huoidities in the lowest layer (nearly saturated conditions 
poleward of latitude 50 degrees) and considerable cloudiness (sore 
than 70 percent total cloud cover) in the polar regions, with a 
secondary aaxiiDUBi of cellular cloudiness at the Equator, but little 
cloud cover in the tropics and subtropics. Surface easterlies appear 
at all latitudes from 70N to 47S, except at the Equator, with north- 
east and southeast trade winds centered at 20N and 27S, respectively, 
and an ITC at 4N, where a mean upward vertical velocity of 4.5 mm s'^ 
closely resembles that found in the first experiment. The westerlies 
aloft show maxima at 31N and 31S (28 and 22 ms"^, respectively) 

at level 6 of the model (nominal pressure, 168 mb). The sea-level 
pressure field displays an equatorial low of less than 1002 mb with bands 
of high pressure (1024 and 1019 mb) at high latitudes (78N and 55S) 
in both hemispheres. In middle latitudes, especially in the Northern 
Hemisphere, both the sea-level pressure field (Fig. 14) and the surface 
winds exhibit large amplitude waves in the first month. 

Within the next few months, the mean January states generated 
by the model in the water planet spin up computation rapidly approach 
the 13-January mean of the first water planet experiment. The evolution 
of the surface pressure profile, for example, as illustrated in Fig. 15, 
shows the change over the first 3 Januaries (see numbers on curves) 
from an initially constant value of 1010 mb to the familiar planetary 
pressure profile, including the equatorial low, subtropical high cells, 
and subpolar lows in both hemispheres. 

When averaged over the last 13 Januaries, the meridional profiles 
of Pg, 0, Vj^, and R for the spin up experiment (not shown) are 
virtually identical to those shown in Figs, 1-4. The corresponding mean 
vertical profiles of U (e.g.. Fig. 5) and T are also nearly identical 
in the two climatologies, so that graphs of the results from the two 
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witer planet experiments plotted or the sane diagram would almost 
coincide; The 13- January mean lynoptic patterns from the spin up 
experiment (see e.g., the se a- k,vel pressure field shown in Fig, 16) 
are also almost indistinguishable from the corresponding fields from 
the first water planet run. 

There are, of course, certain differences between the two water 
planet climatologies which can be seen hy comparing Figs. 12 and 13 
with Figs. 17 and 18, which show the corresponding 15th January 
meridional cross-sections of T and U for the spin up run. The two 
temperature cross-sections are almost identical, but there are some 
differences between the two wind patterns. Despite these small residual 
differences, it does appear from the two water planet experiments that 
the model -generated climate is virtually independent of the initial 
state, Furthermore, the adjustment time in the spin up experiment was 
found to be so short that a horizontally unifoxm state of rest was 
adopted as the initial condition for the remaining perpetual January 
runs. 

Flat continents experiment (002) 

For the next phase of the experiment, geographically realistic 
but topographically flat continents, with zero water storage capacity, 
uniform roughness length (0.3 m) and uniform albedo (0.14), were 
superimposed on the water planet. 

A troublesome problem with all perpetual January experiments 
is the accumulation of snow in the Northern Hemisphere, both on sea 
ice and on land, wherever the Surface temperature is less than 0° C. 

In the January water planet computation, snow did, in fact, accumulate 
from month to month, but only over the sea ice in the Northern Hemi- 
sphere. The snow line in that case was constrained from advancing 
southward by the fact that the sea ice limit in all these experiments 
is arbitrarily fixed. However, no such constraint exists in the model 
on the continental snow line, and, without an annual cycle to provide 
for the melting of snow, it was anticipated that continental glaciation 
over the Northern Hemisphere might spread toward the Equator. A 
15-January run with the flat continent model did indeed indicate such 
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ft glacial advance, as illustrated in Fig. 19, where the snow line is 
seftn to have penetrated as far south as latitude SON over Asia and 
North America after 15 Bonths. However, when the aodel was run for 
an additional 10 Januaries, the snow line not only stabilized but 
even retreated, at least over North /ioerica, as shown in Fig. 20 for 
the 20-th Bonth. In view of the apparently non- catastrophic behavior 
of the continental snow cover, it was considered reasonable to average 
the last 18 months of the 20-January run to derive "climatological" 
maps (Figs. 21-23) for the flat continent experiment, (The reason 
for the apparent constraint on the sovthward spread of the snow cover 
is not immediately obvious, but it appears to be associated with in- 
creased warm advection associated with the development of strong thermal 
gradients between the fixed temperatures of the relatively warm oceans 
and the lowered temperatures of the snow-covered continents.) 

The thermal influence of the continents appears clearly in the 
mean temperature map for the 1000-850 mb layer (Fig. 21), while the 
corresponding monsoonal sea-level pressure distribution (low pressure 
over the excessively heated continents of the Southern Hemisphere and 
high pressure over the cold continents in the Northern Hemisphere) 
is evident in Figure 22. The thermal effect of the continents on the 
500 mb circulation can be inferred from Fig. 23. A comparison of 
the 18-January mean 500 mb contour pattern for the flat continent 
experiment with the 13-January mean for the water planet experiment 
(Fig. 11) shows, first of all, that the 500 mb surface is elevated 
over the heated continents of the Southern Hemisphere and lowered over 
high latitudes of the Northern Hemisphere. Otherwise, the distortion 
of the zonal flow by the continents is rather subtle, the most obvious 
effect being an east coastal trough over North America with a ridge 
in the North Atlantic. 
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Compared with the observed c.Umatological_.January sea-level pressure field 
(Fig, 24), that generated by the flat continent model (Fig. 22) exhibits 
certain deficiencies. The subtropical highs, the Antarctic low, and the 
Aleutian low^ are all too weak, as are the associated pressure gradients, 
compared with the actual climatology. At the 500 mb level, the model 
(Fig, 23) exhibits a more zonal pattern than nature (Fig. 25), especially 
in the Northern Hemisphere, with stronger westerlies in the Arctic. Nota- 
bly absent in the flat continent model climatology at 500 mb are the 
lows over the Canadian Archipelago and eastern Siberia. 

(Mean maps for the flat continent experiment were also computed 
for the last 20 Januaries of the 25-month run, i.e. months 6 through 
25, with essentially the same results, as shown in Figs. 21-23.) 

The influence of the continents on the mean meridional cross-sections 
of temperature and zonal wind appears to be relatively small. The cross- 
sections for the flat continent experiment (not shown) are essentially 
similar to those of the water planet experiments (Figs. 12, 13, 17, 18), 
except for the presence of stronger equatorial easterlies. 

The thermal influence of the continents may be further evaluated 
by comparing the global distributions of precipitation for the water 
planet and the flat continent model. Fig. 26 shows the daily precipita- 
tion rate for the water planet in the 10-th January, while Fig. 27 
displays it for the 20-th January of the flat continent run. In the 
Northern Hemisphere, the continents clearly suppress precipitation in 
winter, while the Aleutian and Icelandic lows enhance it. Over the 
heated continents near the Equator and in the Southern Hemisphere 
precipitation is also augmented, while the oceanic precipitation in the 
equatorial belt is generally reduced compared with that of the water 
planet. The precipitation pattern over South America is quite real- 
istic, but too much rain is computed off the coast of Peru. In the 
Indian and Pacific Oceans, the precipitation maxima are displaced south- 

^A deep Aleutian low does appear after 20 Januaries, but it is located 
over the Bering Strait. Large inter-annual pressure fluctuations are 
found in this region over the 25 Januaries simulated. 
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ward relative to the locations of the Maxima on the water planet, indicating 
that the continents are responsible for the southward displacement of the 
intertropical convergence zone in January over the oceans as well as 
over land* The relatively dry condition in the eastern tropical Atlantic 
Ocean is a particularly noteworthy (and realistic) exas^le of the dis- 
tortion of the zonally symmetric precipitation pattern by the continents, 

The shift of the ITC into the Southern Hemisphere is better illus- 
strated by the precipitation rate map for the 25-th January (Fig. 28), 
showing maxima south of the Et^uator in the Pacific and Indian Oceans. 

However, the dry belt on the eastern side of the tropical Pacific is again 
missing in the model simulation. 

The role of mountains (run 003) 

Both the thermal influences of the continents and the dynamic effects 
of mountains are known to be significant factors controlling the structure 
and circulation of the atmosphere. However, the relative importance and 
quantitative magnitude of each as a climate control is still a subject 
of controversy. As noted, for example, by Kasahara and Washington (1971), 
there is considerable evidence that, with regard to the large-scale climate, 
the thermal effect of continentality is far more important than that of 
orography, although some evidence to the contrary also exists. A cwnparison 
of the January climates generated in the "flat continent" and "mountain" 
experiments with the GISS coarse mesh model should provide some addi- 
tional insight into this problem. 

The orographic experiment differs from those preceding it not only 
in the inclusion of coarse mesh continental topography, but also in the 
use of a dry, isothermal (273 K) initial state. In this case, the model 
is allowed not only to spin up from a state of rest, but also to generate 
its own moisture field (by evaporation from the oceans), as well as vertical 
and horizontal temperature distributions. During the first two months of 
this run, the dry initial state results in extremely low radiative surface 




tea^eratures, especially over the polar night region of the Northern 
Henisphere. The cold continental conditions even produce transient 
snow accuaulations in the subtropics. However, after several Januaries, 
the effects of evaporation from the sea, and the heating of the atmo- 
sphere, begin to establish realistic temperatures over mott of the model 
atmosphere, resulting in a snow line almost identical to that of the 
previous (flat continent) experiment, but with lower teiq>eratures and 
snow over the elevated Antarctic plateau, as shown in Fig. 29 for the 
14-th January. 

One source of difficulty in the mountain experiment is Greenland, 
where the elevated glacier generates extremely low temperatures and, hence, 
excessively high sea-level pressures due to the mc^^thod of reducing station 
pressures to sea level. On the other hand, the snow-covered Antarctic 
plateau produces more realistic surface temperatures on that glacier in 
the January mountain experiment than in the flat continent run, as 
shown in Fig. 30 for the 14-th January. 

Because of the dry initial state in the mountain experiment, the 
model was run for 25 Januaries, and the first 5 months were discarded 
as transients before averaging. The resulting 20-January mean maps 
may then be con 5 >ared with the corre' spending mean maps for the flat con- 
tinent run. 

The 20-January (6-25) mean maps of precipitation rate for the 
"flat continent" run (002) and the "mountain" run (003) are show*^ 
in Figs. 31 and 32, respectively. A comparison of the two fields 
indicates only small effects of the coarse mesh mountains on the pre- 
cipitation in the Northern (winter) Hemisphere. Relatively modest in- 
creases are seen in the region of the Appalachian Mountains (south of 
the Great Lakes) and on the west coast of North America. However, 
much larger differences are found over the warm, moist regions of Africa 
and South America, with greater precipitation indicated in the mountain 
run, not only where the terrain would account for it (e.g., in the Andes 
and over the south African highlands), but also where there should be 
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no orographic component (e.g., eastern South America). Indeed^ the 
change in the orientation and intensity of the precipitation pattern 
over South America appears to be one of the major consequerices of 
the addition of mountains. 

The mean surface temperatures for the last 20 Januaries (6-25 ) , 
shown in Figures 33 and 34 for the flat continent and mountain runs, 
respectively, exhibit not only the expected effects of altitude on tem- 
perature (see, e.g., the Atlantic, Greenland, the Tibetan plateau, and 
southern Africa), but also some residual low temperatures apparently 
resulting from the dry initialization of the mountain run (see, e.g., 
northeastern Siberia). 

Mean meridional cross-sections of zonal winds for the last three 
experiments (Figs. 35-37) illustrate some of the effects of the con- 
tinents and the mountains on the zonal circulation of the model atmo- 
sphere. Fig. 35 is for the water planet run. Fig. 36 for the flat 
continent non, and Fig. 37 for the mountain run. The principal dif- 
ferences between the water planet and flat continent cross=sections 
are found in the equatorial easterlies (stronger with continents) 
and in the horizontal wind shear in the subtropics of the Southern 
Hemisphere (stronger with continents). The effect of mountains, as 
illustrated by the difference between Figs. 36 and 37 is clearly 
smaller than that of the continents. 

(As the experiment is still incomplete, no attempt will be made 
at this time to summarize results or to draw premature conclusions.) 
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Figuws 

1. Me§n January surface pres$ure (ab) vs. latitude for 13 aonths 
(3-15) of aodel run 001 (water planet, with tonal aean 
initial conditions) . 

2. Mean January vertically-averaged tonal wind (as~^) vs. latitude 
for aodel run 001. Initial state, dashed; 13-aonth (3-15) 
average, solid. 

3. Mean January aeridional wind coiqjonent (ca s'^) vs. latitude 

in the lowest layer of the cliaate aodel, run 001. Initial staXe, 
dashed; 13-Bopth (3-15) average, solid. 

4. Mean vertically-averaged January vertical velocity (am s"^) 
vs. latitude for 13 aonths (3-15) of aodel run 001. 

5. Vertical profiles of aean January zonal winds (as~^) vs. pressure 
at latitudes 39N and 37S for aodel run 001. Initial state, dashed; 
13-aonth (3-15) average, solid. 

6. Mean January sea-level pressure field (mb) for month No. 3 from 

model run 001. 

7. Mean January sea-level pressure field (mb) for 13 aonths (3-15) 

of model run 001. 

8. Mean January temperature (K) of layer 1000-850 ab for 13 months 

(3-15) of model run 001. 

9. Mean January temperature (K) of layer 700-500 mb for 13 aonths 

(3-15) of model run 001. 

10. Mean January geopotential height (m) of 700 mb level for 13 months 

(3-15) of model nm 001. 

11. Mean January geopotential height (m) of 5G0 mb level for 13 months 

(3.-15) of model run 001. 

12. Mean meridional cross-section of temperature (°C) for 15-th January 
of model run 001. 

13. Mean meridional cross-section of zonal wind (lO"^ ms”^) for 15-th 

January of model run 001. 

14. Mean Januaiy sea-level pressure field (mb) for first month of model 
run 000 (water planet spin-up experiment). 
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Fi£ures (Cont’d.) 

lb. Evolution of nean January surface '^rossure profile (iib) vs. lati- 
tude during first 3 Months of Model run 000. 

16. Mean Januaiy sea-level pressure field (Mb) for 13 Months (3-15) 
of Model run 000. 

17. Mean Meridional cross-section of tenperature (^C) for 15-th January 
of Model run 000. 

18. Mean meridional cross-section of zonal wind (10’^ ms“^) for 15-th 
January of model run 000. 

19. January snow cover after 15 month run with model 002 (flat continent 
spin -up experiment). 

20. January snow cover after 20 month run with model 002. 

21. Mean January tenperature (K) of layer 1000-850 mb for 18 months 
(3-20) of model run 002. 

22. Mean January sea-level pressure (mb) for 18 months (3-20) of 

mode 1 run 002 . 

23. Mean January geopotential height (m) of 500 mb level for 18 months 

(3-20) of model run 002. 

24. Observed climatological mean January sea- level pressure field (mb) 

(from NOAA and NCAR data) . 

25. Observed climatological mean Januaiy 500 mb geopotential height field 
(m) (from NOAA and NCAR data) . 

26. Mean daily precipitation rate (mm day“^) in January computed for 

the 10-th month of the water planet spin-up experiment (000). 

27. Mean daily precipitation rate (mm day"’^) in January computed for 

the 20-th month of the flat continent experiment (002) 

28. Mean daily precipitation rate (mm day"^) in January computed for 

the 25-th month of model run 002. 

29. January snow cover after 14 month run with model 003 (mountain 
experiment) . 

30. Mean January surface tenperature field (°C) for 14-th month of 
model nm 003. 
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Figures (Cont*d.) 


51. Mean daily precipitation rate (mm day~^) in Januai?/ for 20 months 

(6-25) of model run 002 (flat continents). 

32. Mean daily precipitation rate (mm day'^) in January for 20 months 

(6-25) of model run 003 (mountains). 

33. Mean surface temperature field (°C) in January for 20 months (6-25) 
of model run 002 (flat continents). 

34. Mean surface temperature field (®C) in January for 20 months (6-25) 
of model run 003 (mountains). 

35. Mean meridional cross-secticm of zonal wind (lO"^ ms"^) in January 

for 13 months (3-15) of model run 000 (water planet, spin-up). 

36. Mean meridional cross-section of zonal wind (10"^ ms~^) in January 

for 20 months (6-35) of model run 002 (flat continents). 

37. Mean meridional cross-section of zonal wind (lO”^ ms'^) in January 

for 20 months (6-25) of model run 0Q3 (mountains). 


